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Abstract 

We consider two problems arising in the study of the Schrodinger-Newton 
equations. The first is to find their Lie point symmetries. The second, as an 
application of the first, is to investigate an approximate solution corresponding 
to widely separated lumps of probability. The lumps are found to move like 
point particles under a mutual inverse-square law of attraction. 

1 Introduction 

The Schrodinger-Newton equations is the name given by Penrose [H] to the coupled 
system of equations consisting of the Schrodinger equation for a wave-function ip 
moving in a potential 4>, where is obtained by solving the Poisson equation with 
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source p = \ip\ 2 - Specifically, the system in de-dimensionalised form, is 

~^ + U = (1) 
VV = M 2 - (2) 
and it can be thought of as the Schrodinger equation for a particle moving in its own 
gravitational field. 

Penrose was led to consider this system, or at least its stationary solutions, by his 
reflections on quantum mechanics, [Zj,|H!- The system has been independently con- 
sidered several times jH], PO, sometimes under the name of the Schrodinger-Posson 
system, It is known to be derivable from the Einstein equations of General Rel- 
ativity coupled to a complex scalar field as a slow-motion, weak-field approximation 

Dm- 

The system has been studied numerically [I], [3], and analytically [12] ■ It 
has a conserved energy and conserves total probability. The stationary, spherically- 
symmetric solutions are known numerically [5J, and the lowest-energy state has no 
zeroes and decays exponentially rapidly at infinity. There is evidence that 'lumps' of 
probability can attract each other gravitationally and even orbit each other, but that 
only the ground-state is absolutely stable [2]. 

In this article, we consider two problems arising in the study of the Schrodinger- 
Newton equations. 

The first is to find the Lie point symmetries of the system. This is an investiga- 
tion available for any system of differential equations, and a natural step in the process 
of understanding a system. Some symmetries are clear from the outset: the system 
must have the symmetries of Euclidean space, namely rotations and translations. The 
usual (linear) Schrodinger equation has a symmetry of phase transformation and a 
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symmetry under Galilean transformation and so one expects something similar from 
the Schrodinger-Newton equations. However, as we shall see, the details of the appro- 
priately generalised Galilean transformation contain a reflection of the Equivalence 
Principle from General Relativity. Finally there is a scaling symmetry. 

The second problem is to find an approximation to the Schrodinger-Newton equa- 
tions corresponding to widely separated spherically-symmetric lumps of probability. 
For this, we use a version of the geodesic approximation for the motion of monopoles 
[3]. That is, we use the form of the Galilean transformation and the scaling symme- 
try found in the solution of the first problem to approximate the lumps with rescaled 
copies of the ground-state set into motion, and calculate the conserved energy as 
a functional of the positions, velocities and individual 'probabilities' of the lumps, 
which play the role of masses. Then we treat this energy by the methods of classical 
mechanics. In this approximation, the lumps emerge as point particles moving un- 
der their mutual inverse-square attraction. This can be expected to remain a good 
approximation while the lumps remain widely-separated but one knows from [2] that 
lumps can merge, with the radiation of probability to infinity. 

The plan of the paper is as follows. In the Section 2, we review the method 
of finding Lie point symmetries, closely following ^T] and apply the method to the 
Schrodinger-Newton equations. The details of the calculation are relegated to an 
Appendix, but we present the resulting transformations and interpret them in this 
section. In Section 3, we use the symmetries from Section 2 to write down the solution 
corresponding to a rescaled lump in motion and integrate out the space coordinates in 
the energy. This leads to the classical-mechanical Lagrangian governing the motion 
of widely separated lumps. 
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2 Lie Point Symmetries 
2.1 General Theory 

We choose notation to follow [TT], so suppose we wish to study a system of second- 
order partial differential equations, with independent variables (x l ;i = 1,...,M) 
and dependent variables (u 13 ; (3 = 1, . . . , N). We clump dependent and independent 
variables together as y a and label the equations of the system by an index A: 

H A (x\uP,u%u% m ) :=H A (y a ) = 0, 

where partial derivatives are indicated by subscripts. 

The configuration space is the space coordinatised by y a and a point symmetry is a 
coordinate transformation y~ a = y a (y b ) of the configuration space with H A {y a ) = 0. A 
Lie point symmetry is a point symmetry which arises from an infinitesimal generator, 
a vector field on the configuration space. 

Suppose the Lie point symmetry is written in terms of its infinitesimal generator 



as: 



x n = x n + eXx i + 0(e 2 )=x n + eC(x\u (3 ) + ... 

u a = u a + eXuP + 0{e 2 ) = u a + er ] a {x i ,uP) + ... (3) 



where 



du a 



e=0 d£ 



(4) 

£=0 



de 

then the infinitesimal generator is the vector field X on the space of y a given by 



d d 

X = ^ n 7 — + 7]°' 



dx n du c 
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For a system of second-order PDEs we need to extend the definition of X to the 

second jet-bundle of the configuration space: 

Odd d 
X = — + n a — - + n a +n a — - — 
* dx n ' du a ln du a „ nm du\J 



where 



du, n 



du 



,nm 



i I nm 



e=0 



5c 



e=0 



We call this the prolongation of X but denote it by the same symbol. 
To calculate rj a n and r] a nm we differentiate (J3J): 



(5) 



du c 



du a + edrf + ... 



du c 
dx i 



+ e 



dr] a dr] a dvP 
dx 1 du 13 dx i 



dx l + ... 



(6) 



dx r 



dx n + edC + - 



dx 1 + 



dx 1 du 13 dx r 

and, following [TT], deduce from © using (J3J), (JHJ) and (JJJ) that 



V ,nm + V ,nf3 U ^,m + V ,m/3 U ^n ~ £ ,nm U ,i 



— f k nfi iH ii a 4- n a ii@ — ?/ a — ; 

s ,/37 u ' ,n u ,m u ,k ~ 'I ,/3 u ,nm S l n u ,mi c 

— £ fc /3 (" u "jfc U ,nm + ^n^mk + u<X ,nk U ,m) • 



(7) 



(9) 



The condition that X generate a point symmetry is that the prolonged X satisfy 



XH A = 0, mod H A . 



(10) 



The method will be to solve (fTUj) systematically term by term. 



2.2 Symmetries of the time-dependent Schrodinger-Newton 
equations 

We label variables as 

^yiO -j- /-y» 1 /-y> <"y»2 q , /-y»3 

U 1 = U = 1p, U 2 = V = 1p, U 3 = W = 0, 

and then the time-dependent Schrodinger-Newton equations (JTJ), (J2J) are given by 



H\ := iu )Q + u t n + U j2 2 + ^,33 — uw = (11) 

#2 := ^,0-^,11-^,22-^,33+^ = (12) 

#3 : = U>, 11 + 1W 22 + W33 — 1iU = 0. (13) 

Equation (fTUj) gives the symmetry conditions as: 

XH! = iv u + V u u + V U 2 2 + V u zs-V u ™-V w u = Q, (14) 

Xtf 2 = ivl-Vn-V% i -V v m + V v y> + V w v = 0, (15) 

Xtf 3 = ^+^ + ^-^-^ = 0. (16) 

where these are to hold modulo Ha- 



The method is as follows: For each of the Vj a nm in the system (jl4H16j) . we com- 
pute just the terms that contain those second derivatives u a nm that do not occur in 
(|11 |) -([13 |) . These must vanish separately, which gives rise to the first set of condi- 
tions. Using these conditions, we next calculate all terms in T) a mn containing second 
derivatives of u a , to find a second set of conditions. Finally, assisted by all previous 
conditions, we write the symmetry conditions including all terms, and obtain a final 
set of conditions. This is a straightforward but intricate process which we relegate to 
an appendix, just giving the solution here. 
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There are ten linearly independent point symmetries, including two infinite- 
parameter families. We may list the generators as follows: 

n d d d d n d „ d „ d 
Xi = 2t— + x— + y— + z— - 2u— - 2v— - 2w— , 

at ox oy oz au ov aw 

d d 
X2 = y* — 

ox ay 











X 3 = z a x ^~-- 

ox oz 



d d 

X 4 = z 7> 

oy oz 



d d d d 

X5 = 777, X 6 = — , X 7 = — , X 8 = — , 

at ox ay oz 



x 9 = m)\u^-v^ v 



d_ 

dw ' 



X 
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/ x 1 i 
ajit)—- + -x J 



a'At)(u 



d_ 

du 



^ 1 + ''(*) ^ 
dv J 3 dw 



(17) 



Here a = (a,i(t)) is an arbitrary 3- vector function of time and Q is an arbitrary scalar 
function of time. 

For the first eight transformations, the interpretation is clear. First, Xi is a 
scaling symmetry whose finite form can be written 

t = \- 2 t, x j = A"V, u a = Aw a , 



Under this transformation, the normalisation integral / = J \ip\ 2 d 3 x changes: / = XI. 
We exploit this transformation in the next section to arrive at ()18|1 and Ijl9j) . 
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Next, X 2 , . . . , X 4 are space rotations, and X 5 , . . . , X 8 are translations, in space 
and in time. 

Next, X 9 corresponds to a time-dependent phase-change to ip , with phase-factor 
fl(t), accompanied by the subtraction of Q' from <j). 

Finally, Xio is a time-dependent translation in position, accompanied by a space 
and time dependent phase-change in ip and an additive change to <p. The infinitesimal 
form translates as 

5r = a, Sip = iip(-v.r), 5(f) = — -A.r. 

where v = % and A = The finite form is readily found and we shall need it to 
derive (j!8)) and (^HJ). If A is zero, this is the transformation to a frame with constant 
velocity, in other words a Galilean transformation. The usual Schrodinger equation is 
known to be invariant under Galilean transformation but this calculation shows that 
that property extends to the nonlinear Schrodinger-Newton system. 

If A ^ then we still have an invariance provided we transform the gravitational 
potential. This can be interpreted as a relic of the Equivalence Principle: transforming 
to an accelerating frame is equivalent to adding a constant gravitational field. The 
Schrodinger-Newton equations thus satisfy the Equivalence Principle and we can think 
of this as being due to the fact that they can be obtained by a slow-motion, weak- 
field approximation from the Einstein equations of General Relativity coupled to a 
complex scalar field, jlQj. 

In the next section, we shall use Xi and Xio to put the stationary ground-state 
of total probability one into motion and reduce its total probability. 
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3 The Geodesic Approximation 

We consider an approximate solution of the Schrodinger-Newton equations corre- 
sponding to a number of widely separated 'lumps' of probability, where each lump is 
described approximately by a suitably rescaled copy of the lowest energy spherically- 
symmetric state. (It isn't known, but is rather likely, that the lowest energy spherically- 
symmetric state is the actual ground state.) The lumps are rescaled copies so that 
the total probability is still one. 

From the previous section, we know that each lump, in the absence of all the 
others, can be put into uniform motion by a Galilean transformation. As a station- 
ary solution of the Schrodinger-Newton equations, it is then characterised by its total 
probability, position and velocity. Suppose we have initial data with n lumps with 
corresponding parameters (m^a^, Vj) where Vj = a^ and the overdot stands for d/dt. 
The idea of the geodesic approximation is that, provided the lumps are widely sepa- 
rated, an approximate solution of the Schrodinger-Newton equations can be found by 
assuming that for each lump only these parameters change. To find the equations of 
motion for the lumps, one substitutes the ansatz of several lumps into the energy for 
the theory and integrates out the spatial coordinates, to leave a classical mechanical 
Hamiltonian as a function of (m^ a,, Vj). This is then treated as in classical mechanics 
to give the equations of motion. 

Suppose that ^ (r)e~ lEot is the lowest energy spherically-symmetric state of total 
probability one, taken without loss of generality to be real, and 0o( r ) is the corre- 
sponding potential. Then the rescaled versions, based at a and with velocity v and 



9 



total probability m, are given by: 

1 If 

\P(r, i; m, a, v) = m 2 ^ (m\r - a|) exp(-i(m 2 £ , t - -v.r + - / |v| 2 dt)) (18) 
$(r, t; m, a, v) = m 2 o (m|r — a|). (19) 

For the geodesic approximation, we assume that we have n widely separated lumps 
so that the wave function is 

n 

*(r,f) = X> 

i=i 

n 

= J^r^m^a*,^). (20) 

i=i 

where, to impose the normalisation condition, we require 

n 

E 



i=l 



To find the corresponding potential, we should solve the Poisson equation with 
source equal to l^l 2 , which is in turn a sum. However, since the ground state ipo(r) 
decays exponentially rapidly at infinity and we are assuming that the lumps are widely 
separated, the cross-terms in this sum are everywhere small compared to the diagonal 
terms. Therefore we may take for $ just the sum of the individual potentials due to 
the individual Thus: 



$(r,0 = 

i=i 

n 

= ^$(r,t;mi,ai,Vi). (21) 

i=i 

The energy 8 for the Schrodinger-Newton equations is as follows: 

8 = j QvW.V* + d 3 x (22) 
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with the understanding that $ is obtained from |\1/| 2 via the Poisson equation. 
Now we must substitute \I/ as in (|20p and $ as in (|2Tjl into (|22|) and perform the space 
integrations. To help with this, we need some facts about ipo, 0o and Eq. We have 

\V^ \ 2 d 3 x = —Eq, (23) 
0o|Vo| 2 d 3 x = ^£ , (24) 

|2j3 



|^ | 2 d 3 x = 1, (25) 

where the first two, given in ^2], can be thought of as the Virial Theorem. We also 
need the asymptotic form of 0o( r ) which follows easily from the Poisson equation: 

Mr) = ~ TT1 + h-o. (26) 
47r|r| 

To evaluate (J22)) . we split it into two terms: S = I\ + I2 where 

h = y^V^.V^d 3 ^ 

J 2 = / -§\5Jtyd 3 x 
J 4 

For Ii, we again may neglect cross terms and retain only the diagonal terms, of which 
a typical one is 

y^Wi.V^id 3 x. 

Substituting from (fl"8|) this is 

1 f 1 

—mi I |miV"0o(wi |r — a x |) + -iv 1 ?/' (m 1 |r — a x |) | 2 d 3 r. 
We change the variable of integration to x = mi(r — a x ) and use ()23|) and (|25|) to 



evaluate this as 



1 1 , ,2 

--m^o + g7Tii|vi| . 
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The integral I2 is composed of terms like \ J <&i\tyj\ 2 dx (neglecting terms with prod- 
ucts like ^i^j for i 7^ j because of the exponential fall-off). When % = j, we may use 
(j2H) to find that this is ^mfE . When i 7^ j, we use ()26|) . As long as the lumps are 
widely separated, we may regard the j-th lump as a point particle in the potential of 
the i-th to calculate this term as 

1 mirrij 
167r [a.f — a.j I 

Combining the terms for Ii and I2 and dropping constant multiples of E Q we obtain 

8 ^ 11 16tt V a i - ai . 
i=i i^j ■>' 

Multiplying by 4 and switching from the Hamiltonian to the Lagrangian, we find that 

widely separated lumps move according to the classical mechanical Lagrangian 

1 n 1 

L = -\ m, v i 2 + - > 1 —t- 

2 4:71 ^—^ a; - a; . 
1=1 l^j ■>' 

This is the Lagrangian for point particles of masses m; at positions aj moving un- 
der their mutual, inverse-square law, gravitational attractions: the widely separated 
lumps of probability attract each other with 'gravitational constant' |k 
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Appendix 

Following the strategy of [TT] outlined above, we first write out the parts of i] a mn which 
appear in ([T3|) - (fTT)|) . These are 

V nn = ~ 2£ n U * ni — 2£ pU n U , nk , 
V a n = -2(^0 + e 2 l<12 + ft<13 

+ e>,i<io + e>,i<i 2 + e, u u,iu% 

+ c°, v v,iu% + £>,i< 2 + e, v v,m% 

+ e>.1^10 + + ^^1<13), 

V%2 = -2(^V,20 + £>,21 + £ 3 2 «23 

+ 2« 20 + 2« 21 + £ 3 U M,2«,23 

+ £°^,2M,20 + i\ v V,2U,2l + £V,2« 23 

+ 6°™W2«,20 + £^,10,211,21 + 2M,23), 

^33 = -2(C 3«,30 + ^3 M ,31 + ft" 32 

+ £°« M ,3« 30 + 6^,3^,31 + ^ U M,3M,32 

+ ^^,3^,30 + ^ V'3 U .31 + ^^,3^,32 

+ 3«30 + 3«,31 + £ 2 W W 3 M 32). 

The symmetry conditions (|14|) - (|T^j) . retaining just second order terms, are r] a u +^22 + 
^33 + ••• = 0) which gives the following first set of conditions on £: 

?j> = 0, (27) 

e,j = 0, (28) 
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£,2 ~~ £,1) £,3 — 6^1) £,3 — 6% (29) 

where i = 0, 1,2,3 and j = 1,2,3 (in future, we use i and j as subscripts summed 
over just these ranges). 

Now we shall return to (jHJ) and calculate all second order derivatives contributing 
to T) a n n , using the above first set of conditions. For a = u we find 

Vn = V u ,uU,ii + V u , v v,ii+V u , w w, 11 -2C j il u lj , 

V U 22 = r] U >u U 22 + T] U !V V t22 + T) U >W W 7 22 - 2^U^, 
Via = ^>,33 + ?7>,33 + V U W W,33 ~ 2£\ 3 U 3j . 

Inserting these into (JUJ), considering only the contribution from second order deriva- 
tives, we find 

X#! = Vn + V a 22 + V% + - 

= (Vl + Vl + Vl) 11 + M ,22 + U 33) 
"2 (^,11 + £>,22 + £> 5 3 3 ) + ••■ 

This must vanish mod Hi, which gives us the condition 

c\. = e,2 = ev (30) 

Using this and (J2llj) we obtain: 

£*ll = — ^22 = — ^33' 

s-2 ^2 _ C*2 

~~ S ,11 — S ,22 ~~ S ,33' 

— £^11 = — £ 3 22 = £ 3 33- (31) 
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These ((|30p and ([31))) are the second set of conditions. No extra conditions arise at 
this stage from considering Xi^2 or as they are similar at second order. 

Now we write out the symmetry conditions fully, with all terms, and apply the 
conditions so far obtained. The first is: 

= i(v u ,o + V u ,u u ,o + V w ,v v ,o + V u , w ™,o 

+ <ii + V u 22 + V u 33 ~ V u w - r] w u 

+ 2 (v u ,i u u,i + 7/ V,i + V"iv,w,i + V\u t2 

+ V U ,2vV,2 + V U ,2u,W,2 + < 3U « 3 + rfavj 
+ V U ,3w W ,3) - - ^UU,2 - £ 3 ll« 3 

- ~ & u ,2 ~ f,22 u ,3 - £.33" 1 

- &«,2 - £?33«,3 + V U ,uu («? + M J + M j) 

+ »7% W + « 2 + « j) + ^ W + «> j + Wj) 

+ + U ,2^,2 + «,3f ,3) 

+ *7%, ( U A W ,1 + «,2W,2 + W 3W, 3) 

+ »7%, (",1^1 + V )2 W 2 + ^3^,3) 

+ + + V U ,w W ,H 

+ V U ,u U ,22 + V U ,v V ,22 + V U ,w W ,22 

+ V U ,u U ,33 + V U ,v V ,33 + V U ,w W ,33 

-2e>,u-2e> ,22-2e> ,33 
= 0. (32) 

Equation (|3*2|) must be satisfied identically in the derivatives of w Q mod Ha, so some 
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conditions can be quickly read off. From the coefficients of vp and w^. 



= v u >w = o, 

from quadratic terms in u a f 
from the coefficient of u t j\ 

-<^ + 2^ JU - e, u - e,22 - & = o, 

for j = 1,2,3. The remaining terms in this equation are 

**.o (tu - ft) + «,ii - 2 ft) 
22 « u - 2ft) + n, 33 « u - 2ft) 

+<ll+<22+<33-^-<« + <0 = 0. 

again, mod H\. Using H\ = 0, the terms with 77"^ can be simplified to obt 
iftu Q + 2ftu u + 2ftw 22 + 2ftw 33 - r] u >u uw 

+^ + ^-<ll-<22-<33-<0 = 0, 



mod Hi. Hence 

ft = 2ft = 2ft = 2ft, 

leaving 

(ft - V u , u )uw + rfw + r] w u - V 2 rf - i< = °- 
Conditions (|3^j) and (J53)l give 

rf = uU (t, x, y,z) + U (t, x, y, z). 
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Now, by (|2H1) £° is a function of t only, so by (J35j) the £ J are only first order in 

■3 ■ 



e = 2e(t), 

f 1 = xe'(t)+X(t,y,z), 

e 3 = ze'(t)+Z(t,x,y), (39) 



and (|33J) reduces to 



'C:„ = 277--^. (40) 
Note that £ J at this stage satisfies the conformal Killing vector equation: 

^. + ^. = 2e'(t)5 ij :/../ /(). 

To solve, use conditions (f2~H|) and (J3T|) in (jHOj) and integrate to obtain: 

f 1 = e / (t)x + a 2 (t)y + a 3 (t)z + a 6 + X(t), 

£ 2 = -a 2 (t)x + e\t)y + a 4 (t)z + a 7 + Y(t), 

£ 3 = -a 3 (t)x-a 4 (t)l/ + e'(t)2 + a 8 + Z(t), (41) 

where the functions of time and the constants are as yet unknown, but real. 

Now we consider the other two symmetry conditions, taking into account condi- 
tions on £ we have found above. The same process with Xi/ 2 will give 

it\ = -2t^.„, (42) 
(V v ,v ~ t°o>w + VV - irfp - 7] v w - rf»v = 0, (43) 
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rf = v V(t, x, y, z) + V(t, x, y, z). 
And for XH 3 , in which time derivatives don't appear, 



(44) 



Xtf 3 = ^1+^2 + ^33-^-^ 

= <ii + V% + V w 33 ~ rfv - rfu 

+ V^wj) + r/>,n + 77>,u + 0,n 

+ ?7>,22 + V W V V,22 + V^,22 + 7/>,33 
+ T] W V V >33 + T] W W W >33 - 2£>, n - 2£>, 22 

3 3 W ,33 



which must vanish mod if 3. This gives 

v w ww = v w u = v w v = 0, 
v w 3W = 0, 



- £°o) uv + V 2 rf - rfv - V v u = 0, (45) 



so that 



rf = wW{t) + W{t,x,y,z). 

We will now insert this, (jHHjl and (jUj) into equations (fTTj) . pHj) and (|43|) to calculate 
the rf and £\ Starting with (p?Tj) : 

mu(2e'(*) + W) + u(W - V 2 t/ - iU fi ) +wU- V 2 U - iU = 0, 
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whence 

t/ = o, W = -2e'(t), W = V 2 U + iU fi . (46) 

Next, from (g3J) 

-vw(W + 2e\t)) + v(V 2 V - iV fl -W)-wV + V 2 V + iV = 0, 

whence 

V = 0, HP = VV - iV; . (47) 

And finally 

uv(W - 2e'(t)) + V 2 W - uvU - uvV = 0, 

giving 

U + V = -2e'(t), V 2 W = 0. (48) 

We now have many conditions and there are different routes possible. From (jHH) , (j4*T|) 
and (|IHJ) we obtain 

a 2 = a 3 = fl 4 = 

[/ = l e "(x 2 + y 2 + 2 2 ) + -(xJSf' + yF' + zZ') - 2e + zfi(t) 

where Q is a real function. Now recall that V — U. From ()46)1 we obtain an expression 
for W, while from (|47|) we obtain a different one unless e" = 0. Thus 

1 

e(t) = ait + -as 

in terms of constants a\ and 05. We now have enough to assemble the generators as 



= 2ait + a$, 

= aix + a 2 y + a 3 z + a 6 + X(t), 
20 



e = 
e = 



—a 2 x + aiy + 04,2 + a-j + Y(t), 
—a 3 x — a^y + a\Z + as + Z(t), 



u 



—v 



ir.a' - 2ai + iVt{t) 
-ir.a! + 1a\ + i^{t) 



rf" = -2a x w- -r.a" -fi'(t). 

where a(t) = Z(£)) is an arbitary vector function of time and f2(i) is an 

arbitary scalar function of time. In terms of these the symmetry generator is 



d d 
X = [2ait + a 5 ] — + [axx + a 2 y + a 3 z + a 6 + X (t)] — 

d 

+ [—a 2 x + axy + a^z + a 7 + Y(t)\ 7^ 

d 

+ [— a 3 x - a 4 y + a x z + a 8 + Z(t)] — 

oz 

1 

+u -21". a' 



2ax + iSl(t) 



-ir.a' + 2ax + iCl(t) 
-2a x w - ^r.a" - Q'(t) 



d_ 

du 
d_ 

dv 



d_ 

dw 



which leads to (fTTj) . 
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